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The effect of zinc substitution on the iocai eiectronic structure of several cuprates is investigated 
using first-principles cluster calculations. Clusters comprising 5, 9, and 13 copper atoms in the 
cuprate plane of La2Cu04, YBa2Cu307, and YBa2Cu40s are used. Spin polarized calculations 
with different multiplicities in the framework of density functional theory enable a detailed study of 
the changes in the charge and spin density distribution induced by Zn substitution. Furthermore, 
doping with charge carriers in the above materials is simulated and the resulting changes in the 
charge distribution are compared to the changes induced by Zn impurities. These differences are 
then discussed in terms of a phenomenological model related to properties expected from the generic 
phase diagram. The effects of zinc substitution are rather local and as expected the absolute values 
of the Mulliken charges at both nearest and next nearest neighbor oxygens to Zn are larger than 
in the unsubstituted clusters. The calculated electric field gradient at Cu sites that are nearest 
neighbor to Zn is found to be somewhat larger than in the unsubstituted cluster whereas that of 
next nearest neighbors is about 5 % smaller. We conclude that the satellite peak in the Cu NQR 
spectrum occurring upon Zn substitution in YBa2Cu307 and YBa2Cu408 has its origin at Cu that 
are next nearest neighbors to Zn. 

PACS numbers: 74.72.Bk, 74.72. Dn, 74.25.Jb, 74.62. Dh, 74.25.Nf 



I. INTRODUCTION 

It is well known that the addition of a small amount 
of magnetic impurities destroys superconductivity in 
conventional superconductors but substitution of non- 
magnetic impurities is in general harmless. In hole doped 
cuprate superconductors, however, non-magnetic impu- 
rities like Zn and others are detrimental^ and supercon- 
ductivity disappears with a few percent of Zn impurities. 
This fact has initiated an enormous activity in exper- 
imental and theoretical investigations on the effects of 
impurities that substitute copper in the Cu02 planes. 
In spite of the great experimental and theoretical effort 
which has been done so far, the mechanism of the de- 
struction of superconductivity by Zn impurities in the 
cuprates is still controversial. 

In contrast to the vast amount of theoretical investiga- 
tions on the change in all possible collective phenomena 
that may occur upon Zn substitution, only few informa- 
tions are available about the changes in the local elec- 
tronic structure. Large supercells would have to be used 
in traditional band structure calculations to cope with 
the problems of broken lattice periodicity. Since already 
in the unsubstituted materials, the unit cells contain a 
large number of atoms, such supercell calculations are 
prohibitively complex. 

In the present work, the electronic structure around a 
Zn atom substituting a Cu is studied by ab-initio cluster 
calculations. These methods which are not based on pe- 
riodic lattice structures, are particularly suited to calcu- 
late the changes in the charge and spin density distribu- 
tions that are occurring upon the replacement of a regular 
atom by a new species. Cluster calculations for Zn and Ni 
substituted YBa 2 Cu307 have previously been reported 
by Kaplan et aim using clusters with five copper atoms 



in the Cu02 plane for the pure material. They employed 
M0llcr-Plcssct perturbation theory to investigate on the 
differences in the atomic charges between unsubstituted 
and impurity substituted clusters. Here we present re- 
sults obtained on larger clusters comprising 5, 9, and 13 
copper atoms in the Cu02 plane. Spin-polarized calcula- 
tions within the framework of density functional theory 
are carried out and their accuracy is tested by comparing 
theoretical values for electric field gradients (EFGs) with 
those derived from experiments^. 

It is well known that the doping of the cuprates by 
holes destroys the long-range antiferromagnetic correla- 
tions. On the other hand, it has been proposed that 
Zn substitution enhances local antiferromagnetic correla- 
tions^. Abrikoso\£ referring to the calculations of Ka- 
plan et alwk proposed that Zn, creating an excess positive 
charge at the copper site, reduces the hole concentration 
in its vicinity. He pointed out that this would lead to 
the appearance of a bubble of an insulating spin density 
wave phase around the impurity. We therefore investi- 
gate here in detail the redistribution of charge and spin 
densities around Zn atoms in I^CuOj, YBa2Cu30y, 
and YBa2Cu408- In Sec. El the cluster method will be 
described. In Sees. IIIII and IIVI the resulting changes in 
the spin and charge density distribution induced by Zn 
are given. In addition, in Sec. IIVI simulations of doping 
with charge carriers are performed and the changes in the 
charge distribution are compared to the changes induced 
by Zn substitution. These differences are then discussed 
in terms of a phenomenological model related to proper- 
ties expected from the generic phase diagram. Sec. ]V\ is 
devoted to the calculation of EFGs at Cu sites near Zn 
impurities and their comparison to experiments. Sec. I VII 
contains a summary and conclusions. 
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II. THE CLUSTER METHOD AND 
COMPUTATIONAL DETAILS 

A cluster is a careful selection of a contiguous group of 
ions within a solid. The specific choice of the atoms that 
make up a cluster is such that it allows predominantly 
localized properties of a target atom and its vicinity to be 
calculated. A cluster consists of three regions. The target 
atom and at least its nearest neighboring atoms form the 
center of the cluster and the corresponding electrons are 
treated most accurately using first-principles all-electron 
methods. This core region is embedded in a large cloud 
of a few thousand point charges at the respective lattice 
sites imitating the Madelung potential. Point charges at 
the border of the core region are replaced by basis-free 
pseudopotentials to improve the boundary conditions for 
the electrons in the cluster core. These pseudopotentials 
make up the so-called screening region. 

It would be desirable to have clusters that contain as 
many atoms as possible in the core region but there are 
two computational limitations to the cluster size: the 
available computer resources and the convergence of the 
self consistent field procedure. However, these limits have 
been pushed further since our first use of the cluster tech- 
nique (see Ref^) which now allows to use larger clusters 
and hence to observe effects which could not be seen be- 
fore. 
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La2Cu04 

YBa 2 Cu 3 7 

YBa 2 Cu 4 08 


Cug04 2 /Cui 2 La 5 o 51 663 
Cu 9 033/Cu3oY 16 Ba 16 42 573 
Cu 9 3 3/Cu 2 iYi 6 Bai 6 42 573 
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780 
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La2Cu04 

YBa 2 Cu 3 07 

YBa 2 Cu 4 8 


Cui 3 62 /Cui 2 La 7 4 75 971 
Cui 3 049/Cu38Y 24 Ba 24 6 2 841 
Cui 3 049/Cu 25 Y 24 Ba 2 4 6 2 841 


1313 
1144 
1144 



TABLE I: Compilation of the used clusters and some of their 
defining properties: number of atoms (N), number of electrons 
(E), number of basis functions (B). 

In this work three clusters of different size (with 5, 
9, and 13 copper atoms in the core region) are studied 
and their constitutive properties are given in Table |I| 
For clarity these clusters are labeled X/Y where X is 
the chemical formula of the core region and Y the for- 
mula of the ions of the screening region represented by 
pseudopotentials. Fig. ^ depicts the core region of the 
Cui3062/Cui2La74 cluster and Fig. [^Jshows the planar 
atoms in the core regions of the differently sized clusters. 

The atomic positions were chosen according to crys- 
tallographic structure determinations. For La2Cu04 
we used the tetragonal structure with lattice constants 
a = 3.77 A and c = 13.18 A and with a Cu-O(apex) 
distance of 2.40 A as given in Rcf. 8 . For YBa2Cu307 
and YBa2Cu40g the orthorhombic structure data from 
Ref£ and RefiiS, respectively, were used. For the atoms 



in the core region the standard 6-31 1G basis sets were 
employed. The electronic structures in all the clusters 
described above were determined with the Gaussian 03 
quantum chemistry packaged in the framework of den- 
sity functional theory incorporating the exchange func- 
tional proposed by Beck o 12 ! 1 ? together with the corre- 
lation functional of Lee, Yang, and Parui^ (specified by 
the BLYP keyword in the Gaussian 03 program). We 
note that band structure calculations with the local den- 
sity approximation (LDA) give no magnetic solution for 
La2Cu04. In cluster calculations, however, the spin state 
can be chosen and in all cases (irrespective of the func- 
tional (LDA or generalized gradient approximation)) the 
state with the lowest variational energy is found for the 
antiferromagnetic arrangement (see Refsii^i^). 

The results of each calculation were examined with the 
Mullikcn population analysis which gives a description of 
the charge and spin densities in terms of the individual 
atoms and also the constituent orbitals. Other properties 
such as EFGs and hyperfine fields were also recorded at 
each copper center. 

We have shown previously that with the cluster 
method all contributions to the EFG can be determined 
and the resulting values give a direct comparison to NMR 
and NQR experiments with good agreementiSiiLi^. (For 
planar copper atoms with vanishing asymmetry param- 
eter r\ the connection between the z-component of the 
EFG, V zz , and the quadrupole frequency, vq, is given 
by vq — (e/2h)QV zz , where Q is the nuclear quadrupole 
moment.) Moreover, spin-polarized calculations enabled 
the evaluation of magnetic hyperfine interaction energies. 




FIG. 1: Schematic representation of the core region of the 
Cui 3 06 2 /Cui 2 La74 cluster consisting of 13 copper, 36 planar, 
and 26 apical oxygens, all represented by basis functions. 

One of the great advantage of the cluster method com- 
pared with methods using periodic boundary conditions 
is the ability to study the effect of impurities. Band struc- 
ture methods have to include large supercells to cope 
with the breaking of the translational symmetry. In con- 
trast, the only technical complication emerging in the 
cluster method when introducing impurities is a possible 
reduction of the point group symmetry which sometimes 
requires more computational time. 

Since we want to investigate the changes introduced by 
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substituting a copper ion by a zinc ion we first examined 
the influence of lattice relaxation as the sizes of the two 
dications are known to be different. Therefore, we have 
optimized in the Cu5 26/CusLa34 cluster (for La2Cu04) 
the positions of the four oxygens adjacent to the central 
copper atom. We then repeated the same optimization 
procedure in a cluster where we replaced the central Cu 
by a Zn atom. The result is a small increase of the Zn-0 
distance by 0.0069 A corresponding to a tiny reduction 
of the EFG at the Cu adjacent to Zn of less than 0.5 %. 
We therefore neglect in the following possible lattice dis- 
tortions. 




FIG. 2: Planar copper and oxygen atoms of the core re- 
gion (without apical oxygens) for the clusters comprising 5(a), 
9(b), and 13(c) planar copper atoms. 



III. SPIN DENSITY DISTRIBUTION 

We have performed extensive spin-polarized calcula- 
tions using all possible spin multiplicities. The general 
result of these quantum chemical cluster calculations is 
that they well reproduce the antiferromagnetic exchange 
interaction J between NN copper ions with |J| in the 
range of 160 meV. 

Here we illustrate the antiferromagnetic spin arrange- 
ment with the Cui3062/Cui2La74 cluster for La2Cu04. 
With 13 copper atoms we have the possibility of choosing 
a spin multiplicity with an even value ranging from 2 to 
14 (corresponding to a total spin s = (M — l)/2 ranging 
from 1/2 to 13/2). We performed all the corresponding 
calculations and the state with lowest energy is obtained 
by choosing a spin multiplicity of M = 6. In Fig. [3^ 
the spin density along the Cu-0 bonds in the plane is 
displayed for an unsubstituted Cui3 cluster. In the right 
panel (Fig. 0b) the signs of the Mulliken spin densities at 
the planar Cu sites are shown which exhibit an antifer- 
romagnetic spin arrangement, with adjacent copper sites 
having opposite spins. 

Most of the spin density is provided by the singly occu- 
pied molecular orbital which is also highest in energy. It 
is a linear combination of 3d x 2_ y 2 atomic orbitals (AOs) 
on the coppers (n(3d x 2_ y 2) = 79%) and 2p CT AOs on the 
oxygens (n(2p CT ) = 21%). The square of the 3d a; 2_ y 2 AO 
has maxima at distances of 0.345 A from the nucleus but 
vanishes at the nucleus. In general, the spin density at 
a nuclear site is due to s AOs giving rise to the Fermi 
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FIG. 3: (a) Spin density along the 5 Cu and 6 O in the 
Cui3 062/Cui2La74 Cu cluster of spin multiplicity M=6; (c) 
spin density with Zn substitution and multiplicity M=5. 

contact term. In particular, the spin densities at the Cu 
nuclei (cusps in Fig. |2t) can be correlated with the Mul- 
liken spin densities (whose signs are given in Fig. |3Jd) at 
the same and the nearest neighbored Cu nuclei. This 
observation can be used to split the Fermi contact term 
into on-site and transferred terms as has been discussed 
in detail in Refsii£*it We just note here that the on- 
site term is negative if the Mulliken spin density at the 
Cu atom under consideration is positive and vice versa 
and that the transferred term is positive (negative) if 
the Mulliken spin densities at the neighboring Cu atoms 
is positive (negative). Near the oxygens, the spin den- 
sity is provided by the 2-p a AO and the contact density 
is very small for antiferromagnetically arranged copper 
neighbors. 

The substitution of Zn at the cluster center introduces 
an extra electron into the cluster and it is to be noted that 
the spin state of lowest energy is now the M = 5 state. 
The corresponding spin density along the Cu-0 bonds 
(see Fig. shows that the antiferromagnetic arrange- 
ment virtually remains intact. As expected the Mulliken 
spin density on the Zn is very close to zero (Fig. 01). Note 
that now the contact densities at the copper nuclei adja- 
cent to Zn are reduced in comparison to the same sites 
in Fig. [3^ because there are transferred hyperhnc fields 
from three ions only instead of four. These changes in 
the hyperfine fields will be further discussed in Sec. IV Bl 

IV. CHARGE DENSITY DISTRIBUTION 

The cluster method enables a detailed analysis of 
the local charge distributions and their changes upon 
Zn substitution as well as upon simulations of dop- 
ing. We first present the differences observed between 
the Cui2Zn062/Cui2La74 and the Cui3062/Cui2La74 
clusters for La2Cu04. The most prominent modifica- 
tions concern the oxygen atoms adjacent to Zn. We 
therefore then concentrate on the Mulliken charges of 
these oxygens and their differences in undoped and hole 
doped La2Cu04 and in optimally doped and underdoped 
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YBa2Cu307- These quantitative results are then further 
discussed in terms of a qualitative model of extended 
charges. Finally we discuss the implications of our results 
for models which relate the suppression of T c upon Zn 
substitution with enhanced antiferromagnetism around 
Zn substituted regions. 

We first discuss the differences in the charge distribu- 
tions at all atoms. Fig. 0] illustrates the difference in the 
total electron densities in the C11O2 plane between the Zn 
substituted and the unsubstituted La2Cu04 using calcu- 
lations on the corresponding CU13 clusters (with the cen- 
tral Cu atom being substituted by Zn). In both clusters, 
the multiplicity leading to the lowest ground state energy 
was chosen: M = 6 in the unsubstituted and M = 5 in 
the Zn substituted cluster. The framed area insinuates 
the core region of the CU13 cluster. The solid (black) con- 
tours are positive electron density differences, which are 
more prominent than the dashed (red) contours of nega- 
tive electron density differences since the central Zn atom 
has one electron more than the Cu atom it replaces. The 
white space in the corners of the framed area is signifi- 
cant since the four coppers and 24 oxygens in this region 
are virtually unaffected by the zinc substitution and do 
not show up on the contour map (|Ap| < 0.005 a^ 3 ). 

The difference between the center atoms can be ex- 
plained as follows. As expected the 3d x 2_ y 2 AO of the 
Zn is almost fully occupied which shows up in the pro- 
nounced black density difference in the center of Fig. 0] 
The red contour lines are concentrated on 3d xy orbitals. 
This is due to the fact that the 3d xy -electrons of the Zn 
ion are closer to the nuclei than the ones of the Cu ions 
since for Zn (l/r 3 ) = 9.39 a~ 3 but (l/r 3 ) = 7.94 a" 3 
for Cu. Therefore, the electron densities at the border of 
the Zn ion are smaller than at the border of the Cu ion 
which results in red equidensity lines. 

The central zinc atom of course shows the very large 
positive electron density difference but some has been 
transferred to the adjacent oxygen atoms and also to one 
side of the nearest neighbor copper atoms. Thus there 
are significantly more electrons in the 2p CT AO of the oxy- 
gens adjacent to Zn than in the unsubstituted case. The 
occupancy of the 2p CT orbital changes from 1.67 to 1.74 
and the corresponding Mulliken charges are reduced from 
-1.63 to -1.71. 

For the NN copper atoms we observe an asymmetric 
charge distribution. The total charge on these Cu, how- 
ever, almost equals that in the unsubstituted compound. 
On the four coppers which are next nearest neighbors 
to the Zn, there are more electrons in the 3d x 2_ y 2 AO 
compared to the bulk values. These charge differences 
on the NN and NNN copper atoms imply changes in the 
electric field gradients which will be discussed in detail 
in Sec. IVH 

The calculations in RefpS, yield similar changes in the 
Mulliken charges at the oxygens adjacent to the Zn im- 
purity but the changes at the adjacent copper sites are 
much larger than those calculated in this work. This dif- 
ference might be due to the use of bare point charges 




FIG. 4: (color online) Electron density difference contour 
plots in the cuprate plane between the Cui2Zn062/Cui2La74 
and the Cui3C>62/Cui2La74 clusters on a area of 5 x 5 lattice 
constants. The tick marks correspond to the planar copper 
position. The solid density lines (black) occur when the local 
electron density of the Zn substituted cluster is larger than 
that of the corresponding density of the unsubstituted CU13 
cluster. Dashed lines (red), on the other hand, denote that 
the local electron density is larger in the unsubstituted clus- 
ter. Two neighboring lines (solid and dashed) are separated 
by an electron density difference of 0.01 e/a.%. 



in the immediate vicinity of the border of the cluster. It 
was shown in Ref. 19 that especially positive point charges 
should be replaced by pseudopotentials since they present 
a too strong Coulomb potential for the electrons in the 
core region. 

It is remarkable that the effects of Zn substitution are 
rather local as concerns the electronic structure. This 
local screening, however, tells nothing about collective 
phenomena like scattering where it can produce a much 
more extended range of influence. 

We have also determined the electron density dif- 
ferences for clusters representing the YBa2Cu307 and 
YBa2Cu40g materials. The results are very similar 
to those for La2Cu04 shown in Fig. 0] and discussed 
above. Furthermore, to simulate the effect of hole 
doping in La2Cu04 we have repeated all calculations 
for the corresponding clusters Cui30e2/Cui2La74 and 
Cui2Zn062/Cui2La74 with one electron less. These sim- 
ulations correspond to a hole doping of roughly 8 %. In 
the same way, the effects of underdoping YBa2Cu307 
were studied by repeating all calculations with one elec- 
tron more. To discuss the essential changes we concen- 
trate on the Mulliken charges of the four oxygen atoms 
adjacent to Zn or, in the unsubstituted compounds, Cu. 
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The calculated Mulliken charges of the planar oxygens 
in La 2 Cu04 and YBa 2 Cu 3 07 are —1.64 and —1.57, re- 
spectively, as illustrated by the white rectangles in Fig.03 
Both these values are reduced to —1.71 and —1.64, re- 
spectively (diagonally hatched rectangles), for oxygens 
adjacent to substituted Zn atoms. The simulation of hole 
doping in La2Cu04 results in a Mulliken charge of —1.61 
without Zn (horizontally hatched rectangles) and —1.70 
(vertically hatched rectangles) on the oxygens adjacent 
to Zn. An analogous simulation but with electron doping 
was performed for YBa 2 Cu307. If an electron is added to 
the Cui3049/Cu3sY 2 4Ba24 cluster the Mulliken charges 
on the oxygens reduces slightly to —1.58 in the cluster 
without Zn and is virtually unaffected (—1.64) in the Zn 
substituted cluster. 

These quantitative results for the local changes near 
substituted Zn and for the effects of simulating about 
8 % hole doping are now put in perspective with the ex- 
tended changes expected from the generic phase diagram. 
The following discussion is therefore necessarily qualita- 
tive and completely phenomenological since we do not 
consider band structure. It serves, however, to rational- 
ize the main changes that are locally produced by Zn 
substitution. 

The conventional wisdom in cuprate superconduc- 
tors is that the pure parent compounds La 2 Cu04 and 
YBa 2 Cu3C>6 are insulators and that upon "hole doping" 
the missing electrons come from the planar oxygens. The 
holes thus introduced on the oxygens destroy rapidly the 
long range antiferromagnetic order and the materials be- 
come conducting. Although we have pointed out else- 
wher e) 20 ! 21 that this picture is too naive and may easily 
lead to wrong conclusions, a semiquantitative discussion 
on the values of the oxygen charges and their changes 
with hole and electron doping both in the unsubstituted 
and in the Zn substituted compound is instructive. We 
have therefore indicated in Fig. [S] a schematic separa- 
tion into a metallic region (gray) and an insulating region 
(white) with a border line at the arbitrary chosen value 
of -1.62. 

This puts the Mulliken charge on the oxygen in un- 
doped La 2 CuC>4 (white rectangle in Fig. O into the insu- 
lating region and that obtained with hole doping in the 
more metallic region. From simulations of underdoped 
YBa 2 Cu3C>7, however, we find a small decrease of the 
Mulliken charge which still is the conducting zone. 

In each of the four cases considered, Zn substitution 
drives the Mulliken charges at the oxygens deeper into 
the insulating region. In particular, Zn substitution 
in YBa 2 Cu 3 07 reduces pm{0) to a value (diagonally 
hatched rectangles) which is almost equal to pm(0) in 
undoped and Zn free La 2 Cu04 (white rectangle). Hence, 
we observe that the electronic structure in the immedi- 
ate vicinity of a Zn impurity in hole doped cuprates is 
similar to that found in the Zn free and undoped parent 
compounds. 

This schematic description of the effect of Zn substitu- 
tion and charge carrier doping can be interpreted in terms 
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FIG. 5: Schematic view of the Mulliken charge on the oxy- 
gens in the unsubstituted compounds (white rectangles) and 
in Zn substituted compounds adjacent to a Zn atom (diago- 
nally hatched rectangles) for La2Cu04 and YBa2Cu3C>7. The 
results for clusters doped with one hole for La2Cu04 and one 
electron for YBa2Cua07 are plotted with horizontally (in un- 
substituted clusters) and vertically hatched (in Zn substituted 
clusters) rectangles. An approximate border line at —1.62 is 
drawn to illustrate the separation between the insulating an- 
tiferromagnetic phase and the conducting phase. 



of the generic phase diagram: Zn substitution overcom- 
pensates the hole doping effects in its vicinity and thereby 
creates an environment where the antiferromagnetic cor- 
relations are restored. The prediction, however, whether 
this leads to a static or dynamic spin ordering cannot be 
given with our cluster method. 

Several experiments in Zn substituted cuprates have 
also been interpreted as indications for locally enhanced 
antiferromagnetic correlations (see Refi^) induced by 
Zn. This has been emphasized in the work of Abrikosov 6 
who considered the existence of antiferromagnetic bub- 
bles which surround Zn, Ni, or Li impurities in the Cu0 2 
plane on the basis of the general theory of metal-insulator 
transition in layered cuprates. 



V. SITE ASSIGNMENTS IN COPPER NMR 
AND NQR SPECTRA 

Numerous NMR and NQR experiments on Zn substi- 
tuted cuprates have been reported. In the present work 
we refer in particular to the measurements of Itoh et 
al& and Williams et alA who presented 63 Cu quadrupole 
resonance frequencies and spin-lattice relaxation rates 
in unsubstituted and Zn substituted YBa 2 Cu 3 07 and 
YBa 2 Cu40g. It has been observed that the main cop- 
per line broadens and slightly shifts to higher frequen- 
cies when Zn is substituted and that a new Cu satellite 
resonance occurs at a somewhat lower frequency. The 
intensity of this satellite peak grows with increasing Zn 
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concentration. 

In addition, Itoh et alA found that the Cu spin-lattice 
relaxation time, 63 Ti, of the satellite peak is shorter than 
that of the main signal. This was interpreted as resulting 
from locally enhanced antiferromagnetic correlations in 
the neighborhood of Zn based on the site assignment that 
the NQR satellite peak comes from Cu atoms that are 
nearest neighbors to a substituted Zn atom. 

In this chapter we propose a different site assignment. 
We claim that the satellite peak is due to copper atoms 
that are next nearest neighbors to a substituted Zn atom. 
We will provide evidence for this in Sec. IV Al by a di- 
rect calculation of copper EFGs in clusters represen- 
tative of unsubstituted and Zn substituted La 2 Cu04, 
YBa 2 Cu 3 7 , and YBa 2 Cu 4 8 . In Sec.|VB]we will con- 
sider measurements of copper spin-lattice relaxation rates 
in unsubstituted and Zn substituted cuprates to give ad- 
ditional independent arguments that support the site as- 
signment suggested above. Furthermore, in Sec. IV CI a 
detailed investigation about the sensitivity of the EFGs 
on the charge distribution is given. 



A. Electric Field Gradients 

The extra calculations needed to study the changes of 
the copper EFGs in substances doped with impurities 
have been performed - for efficiency reasons - in clusters 
comprising nine copper atoms in the plane. As already 
stated in Sec. lIVI this is justified since the changes in the 
electronic structure upon Zn substitution occur predomi- 
nantly in the immediate vicinity of the impurity (see also 
Fig.©. 

We have evaluated the electronic structure in clus- 
ters for La 2 Cu04, YBa 2 Cu30y, and YBa 2 Cu40s and 
recorded the copper EFGs at all copper sites in the clus- 
ter which slightly differ due to finite size effects. Next we 
repeated the above calculations for Zn substituted clus- 
ters where the copper in the middle of the cluster (site 
M in the inset of Fig. [fjj) was replaced by a zinc. The 
copper EFGs were again calculated at all sites and were 
compared to those at the same sites in the unsubstituted 
clusters. 

The relative differences between the copper EFGs in 
the Zn substituted and in the unsubstituted clusters are 
shown with thick bars in Fig. |SJ The solid (dotted) thick 
lines refer to situations where the copper in the substi- 
tuted cluster is a nearest neighbor (NN) (next nearest 
neighbor (NNN)) to the Zn impurity. The EFG values 
for the NN coppers are all found to be larger than in the 
unsubstituted compounds whereas those for the NNN are 
smaller. 

For consistency, we repeated the above calculations but 
with Zn impurities in the edge and corner position (site 
E and site C). We included the results of these calcula- 
tions with thin bars in Fig. and - as a guide to the eye 
- grouped the copper EFGs resulting from coppers with 
the same distance (NN or NNN) to the Zn impurity by 
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FIG. 6: Relative differences between the copper EFGs in the 
Zn substituted and in the unsubstituted clusters with solid 
(dotted) lines referring to copper atoms which are NN (NNN) 
to the Zn impurity. In the inset: schematic representation of 
the Zn sites in the CU9 cluster with the site labelings used 
in the text (M: middle, E: edge, C: corner). Thick bars are 
for situations where the Zn is replaced at site M whereas thin 
bars denote results from clusters where the impurity has been 
inserted at sites E or C (see text). 



gray blocks. It is to be expected that asymmetrical sub- 
stitution leads to somewhat deviating values, however, 
the qualitative findings of the above paragraph remain 
unchanged. 

We therefore conclude that the observed satellite peak 
with lower NQR frequency originates from copper nuclei 
that are NNN to Zn. The NN will only broaden the 
main line and induce a shoulder towards higher frequency 
which is clearly visible in the Cu NQR spectrum^. 

Note that in Rcfs.M it was suggested by probabilistic 
arguments that the satellite peaks may originate from 
copper sites with one Zn NN. However, this site assign- 
ment is not unambiguous, since the probability of finding 
a copper coordinated with one Zn NN is the same as that 
of the Zn being in a NNN position to the copper, at least 
for low Zn concentration. 



B. Comparison of the calculated hyperfine field 
with the NMR spin-lattice relaxation rate 

In this paragraph we will give additional independent 
arguments that support the above conclusion about the 
site assignments in 63 Cu NQR spectra. 

The spin-lattice relaxation rate with the applied field 
along the c-axis, 63 T^r, is due to fluctuations of the hy- 
perfine fields in the perpendicular directions, i.e. in the 
Cu0 2 plane. Adopting the notation from the Mila-Rice 
spin Hamiltonian^, there is an on-site hyperfine interac- 
tion at the copper nucleus, denoted by A±, and a trans- 
ferred field B. For the case of fully antiferromagnetically 
correlated copper spins, the rate 63 T 1 ~ 1 , is proportional 
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to 

63 T^} oc (A ± - 4B) 2 (1) 

since each of the four NN contributes one B. A Zn NN, 
however, will hardly transfer any spin density to the Cu 
(see Fig. |3J) and we expect 

63 T^ oc (Ax - 3B) 2 . (2) 

Typical values^ are A±_ = 0.29 j«eV and B = 0.4 ^teV, 
which would reduce ^T^ 1 by 1 - (A x - 3S) 2 /(A ± - 
4S) 2 = 52 %. Note that if there were no correlations at 
all, this reduction would still be 1 - (A\ + 3B 2 )/(A 2 j_ + 
AB 2 ) = 22 %. 

To obtain a more quantitative statement we can 
use the results of our cluster calculations in the 
Cug042/Cui2La5o and CusZnC^/Cui^Laso clusters for 
the La2CuC>4 system. The circles in Fig. [7| indicate the 
calculated Fermi contacts, D(Cu), in the unsubstituted 
Cug cluster recorded at different positions in the clus- 
ter (see inset of Fig. |SJ with different number of nearest 
neighbors. The slope of the fitted straight line gives the 
transferred contribution and the intersect the on-site con- 
tribution. The squares denote results in a Zn substituted 
cluster. The arrows are meant to indicate that upon los- 
ing one Cu neighbor the Fermi contact drops by about 
0.6-B. Since there is some spin density on the oxygen 
between Zn and Cu, the change of the total transferred 
spin density is not from AB to 3B but to SAB which 
still leads to an enormous reduction of the spin-lattice 
relaxation rate. 

For a Cu nucleus adjacent to a Zn we therefore ex- 
pect a considerably reduced spin-lattice relaxation rate. 
If these copper nuclei contributed to the satellite peak 
they would relax slower than those contributing to the 
main peak in contrast to the observations (see Fig. 7 
in Ref£). In addition the spin-lattice relaxation rate of 
nuclei contributing to the main peak would not change 
when substituting Zn atoms in the cuprate plane. How- 
ever, the data in Fig. 6 of Ref£ show that, in fact, the 
relaxation of the copper nuclei in the main peak is slower 
in Zn substituted compounds. This behavior is expected 
only if, instead, the Cu atoms nearest neighbored to Zn 
contribute to the main peak. 

We thus conclude that the satellite peak at which a 
larger rate is observed cannot be attributed to a Cu which 
is nearest neighbor to Zn thus providing an additional 
argument to the assignment made in Sec. IV Al 



C. Sensitivity of EFG values on charge distribution 

A detailed description of the evaluation of EFG in the 
framework of cluster calculations is found in Refpi^. Here 
we just point out the relevant aspects that are connected 
to the rearrangement of charge distributions in the neigh- 
borhood of a Zn impurity. In cluster calculations, the 




number of nearest copper neighbors 



FIG. 7: Contact spin density D(Cu) in a^ 3 vs. the number of 
nearest copper neighbors. The circles result from calculations 
for the unsubstituted compound whereas the squares refer to 
Zn substituted clusters. 

charge distribution is determined by the occupied MOs. 
The m th MO is represented as a linear combination 

n n tlk 

<Ur) = £ Rk) = E X ^' k B K ,k(r- Rk) 

K=\ K=l fc=l 

(3) 

of riK atomic basis functions Bk,u centered at the nuclear 
sites K = 1, . . . , n, and the c^' k are the MO coefficients. 
We assume in the following that the target nucleus Kt 
is at Rk t = 0. The contribution of the MO (f> m to the 
EFG at Kt is given by the matrix element 

^ = (^(01 % \<t>m{f)) 

n n riK Ul 

= X/ X/ E X/ C m' kc m l x 
K=l L=l k=l 1=1 

ix l x^ — r 2 5^ 

x{B K<k (f-R K )\ — I B L .i{r-R L )). (4) 

Thus we can identify three types of contributions: (i) 
on-site terms from basis functions centered at the target 
nucleus (K = L = Kt), (ii) mixed on-site and off-site 
contributions, and (Hi) purely off-site terms with K ^ 
Kt and L ^ Kt (see Ref4& for further details) . 

To rationalize the changes of the EFG on the copper 
sites upon replacing a Cu by a Zn atom, it is sufficient to 
concentrate on the MOs which contain partially occupied 
3d x 2_ y 2 and 3d 32 2_ r 2 AOs at the target nucleus. The 
EFG is then given by the contributions R from all other 
orbitals and these two frontier orbitals: 

V ZZ =R+ ^{N 3z 2_ r 2(r- 3 ) 3z 2_ r 2 - N x 2_ y 2 ( r~ 3 ) x 2_ y 2) 

(5) 
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Here, N are the occupancies of the corresponding or- 
bitals. In Table ITD we collect the values obtained at the 
NNN Cu in the pure and the Zn substituted clusters, 
respectively. We illustrate this fact with the La2Cu04 
substance and note that the contributions from other or- 
bitals (R) is marginally changed upon Zn substitution. 

The decrease of 4 % of V zz at the NNN Cu is thus 
due to an increase in N x 2_ y 2 by 0.010 and a decrease 
in N 3z 2_ r 2 by 0.002. This demonstrates that EFGs are 
rather sensitive to deviations of the electronic charge den- 
sities. 

The fact that the partial occupancies N of the AOs 
are very similar to the partial Mulliken charges pm which 
are also included in Table [H] gives evidence that despite 
their simplistic concept the Mulliken charges can be sig- 
nificant quantities in the analysis of charge distributions 
in cuprates. 

unsubstituted Zn substituted 

3d x 2_y2 3d^ z 2_ r 2 3d x 2_y2 3d$ z 2_ r 2 

N 1.334 1.921 1.344 1.919 

Pm 1.377 1.931 1.387 1.929 
( r~ 3 ) 8.058 8.004 8.048 7.999 

TABLE II: Occupancies N, partial Mulliken charges pM, and 
expectation values { r -3 } for the relevant copper AOs in un- 
substituted (M = 2) and Zn substituted (M = 1) Cug clusters 
representative of La2Cu04. 



VI. SUMMARY AND CONCLUSIONS 

The objective of this work was a first-principles inves- 
tigation on the effect of Zn impurities that substitute Cu 
atoms in the CuC-2 plane of the cuprates on the local elec- 
tronic structure. This was performed in the framework 
of spin-polarized density functional theory with localized 
basis functions in clusters representative of La2Cu04, 
YBa 2 Cu 3 07, and YBa 2 Cu40 8 and the respective Zn sub- 
stituted compounds. As a significant improvement of our 
previous worki^ii we were able to include up to 13 pla- 
nar copper atoms together with the surrounding planar 
and apical oxygen atoms in our clusters. 

For the unsubstituted compounds the antiferromag- 
netic spin arrangement with the Mulliken spin densities 



of neighboring copper atoms having opposite sign could 
well be reproduced. Calculations in clusters where a Cu 
was replaced by a Zn atom revealed almost no spin den- 
sity at the Zn site whereas the antiferromagnetic spin 
alignment around the impurity was left virtually intact. 

Using the same clusters, we observed in general that 
the changes in the charge distribution upon Zn substi- 
tution remain rather local and does not extend to more 
than the next nearest copper atoms. In particular the 
Mulliken charge at the oxygens adjacent to the Zn im- 
purity was significantly reduced to more negative values 
compared to the Zn free case. We thus obtain a reduced 
"hole concentration" in the immediate neighborhood of a 
Zn impurity, i.e. on the four next nearest oxygen atoms, 
i.e. on the four next nearest oxygen atoms. This result 
is in agreement with that obtained previously by Kaplan 
et al£. For a comparison of the effects of Zn substitution 
with those of hole doping, we then simulated a hole dop- 
ing of about 8 % in unsubstituted La2Cu04 by removing 
one electron from the Cui3 cluster and observed that in 
general impurity substitution seems to compensate the 
effects of hole doping on the charge distribution. We in- 
terpreted this as evidence for enhanced antiferromagnetic 
correlations around Zn impurities. 

Furthermore, we also calculated electric field gradients 
at Cu sites in unsubstituted and Zn substituted clusters 
and found that the EFGs at Cu sites which are nearest 
neighbors to a Zn impurity are slightly larger than in the 
unsubstituted compounds whereas those at sites NNN to 
Zn are about 5 % smaller. We therefore proposed a new 
site assignment for Cu NQR spectra in Zn substituted 
compounds: the satellite peak at lower frequencies orig- 
inates from NNN Cu to Zn whereas the NN Cu induce 
a shoulder in the main peak towards higher frequencies. 
This site assignment is corroborated by arguments con- 
cerning the different spin-lattice relaxation rates of nuclei 
in the main and the satellite peak. 
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